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We study the partition of energy between waves and vortices in stratified turbulence, 
with or without rotation, for a variety of parameters, focusing on the behavior of the 
waves and vortices in the inverse cascade of energy towards the large scales. To this end, 
we use direct numerical simulations in a cubic box at a Reynolds number Re » 1000, with 
the ratio between the Brunt-Vaisala frequency N and the inertial frequency / varying 
from 1/4 to 20, together with a purely stratified run. The Froude number, measuring the 
strength of the stratification, varies within the range 0.02 < Fr < 0.32. We find that the 
inverse cascade is dominated by the slow quasi-geostrophic modes. Their energy spectra 
and fluxes exhibit characteristics of an inverse cascade, even though their energy is not 
conserved. Surprisingly, the slow vortices still dominate when the ratio N/f increases, 
also in the stratified case, although less and less so. However, when N/f increases, the 
inverse cascade of the slow modes becomes weaker and weaker, and it vanishes in the 
purely stratified case. We discuss how the disappearance of the inverse cascade of energy 
with increasing N/ f can be interpreted in terms of the waves and vortices, and identify 
three major effects that can explain this transition based on inviscid invariants arguments. 


1. Introduction 

The interaction between vortices and waves is a long-standing problem in fluid turbu¬ 
lence. It can lead to a self-sustaining process that is dominant, for example in pipe flows, 
due to the strong wall-induced shear. A quasi-linear analysis around a mean but unsteady 


flow allows for a prediction of large-scale coherent structures in such flows (McKeon & 


Sharma|2010 

Farrell & Ioannou|2012 

Sharma & McKeon|2013 

I, to the prediction of jets 

in baroclinic turbulence in the atmosphere of planets ( 

Farrell & loannou 

2009 

, and it 


can also be used as a tool to control the onset of turbulence using traveling waves (Moar- 
ref fc Jovanovic|[2010 ). The dynamics of geophysical flows involves complex interactions 
between nonlinear eddies and waves due to the combination of planetary rotation and 
density stratification as well as the formation of shear layers. Typical parameters for 
the atmosphere and the ocean at large scales (larger than about 1000 km) are such 
that the dynamics of the fast waves is slaved to that of the eddies, and this leads to a 
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quasi-geostrophic quasi-linear regime characterized by a balance between pressure gradi¬ 
ent and the Coriolis and gravity forces ( Charne;^|1971 1. Augier & Lindborg ( |2013 ) have 
shown using general circulation models (GCMs) that the ratio between the energy in 
the balanced vortices and the energy in the horizontally divergent motions (including 
internal waves) was on the order of 100. Then, the eddy turn-over time is much larger 
than typical wave frequencies: the Brunt-Vaisala frequency N or the inertial frequency 
/ = 2f2 (with n the rotation rate). For instance, frequency spectra obtained from in- 
situ measurements in the ocean show distinct power-laws for internal waves in the [/, N] 
band and geostrophic eddies at lower frequencies (Ferrari & Wunsch 2009). The quasi- 


geostrophic regime is well understood theoretically and provides a good description of 
the large scales of the atmosphere and the ocean ( Vallis|[^006 e.g.). However, the range 
of scales in these geophysical flows before dissipation prevails at small scales is such that 
other scale-dependent regimes can arise in which turbulence comes into play, with the 
eddy turn-over time becoming comparable to the Brunt-Vaisala frequency, the inertial 
frequency, or their combinations which arise in the dispersion relation of internal waves. 
|UII^ p^83l ), in particular, has underlined the emergence of a regime, coined stratified 
turbulence, between quasi-geostrophic turbulence and homogeneous isotropic turbulence 


(see also the scale analysis by Riley et al. (1981)), and its importance for the mesoscales 
(with horizontal scales between 1 and 100 km, approximately) in the atmosphere. Indeed, 
the measured horizontal spectra of kinetic and potential energy in the upper troposphere 


and lower stratosphere (Nastrom et al. 1984) has a long history of opposite interpreta¬ 
tions: some have claimed that this unambiguous range corresponds to an inverse 

cascade of quasi-2D eddies ( Gage|1979 Lilly|1983 Metais et al. 1996), while others have 
argued that it was in fact a direct energy cascade dominated by internal waves (Dewan 


1979). This example emphasizes the importance of analyzing data — be it observational. 


experimental or numerical — in terms of eddies or vortices on the one hand, and waves 
on the other hand. In this paper, we shall try to disentangle these two types of motion in 
the idealized context of direct numerical simulations of stratified flows with and without 
rotation in a periodic box, focusing on the classical quantities characterizing the statis¬ 
tical behavior of nonlinear interactions in turbulent flows, in particular energy spectra 
and transfer functions. 

One way to distinguish between waves and eddies is through a decomposition of the 


flow onto the normal modes of the linearized system (Bartello 1995 see also below. 
Normal mode decomposition is a standard technique which has been used in a 


& Bartello 2004 

Herbert et al. 2014) and numerical studies (|Bartello 

1995 Waite & 

Bartello|2006a|6 

ISukhatme & Smith|2008 Kitamura & Matsuda||2010 

iurien & Smith 


2012 2014 Brunner-Suzuki et al. 2014), for the Boussinesq equations, which we shall 


investigate here, as well as for other dynamical equations. Recently, another approach 
has been proposed (Biihler et al. 2014), which relies on a few additional assumptions 


in order to be able to achieve a similar decomposition without a complete knowledge of 
the dynamical fields. This is particularly useful when working with observational data. 


such as ship track data or airplane measurements. Gallies et al. (2014) have used this 


method to claim that observational data supported the internal wave interpretation of 
the forward cascade of the atmospheric mesoscales, although Lindborg (2015) argues 


using a similar technique that the forward cascade of energy could in fact be of another 


type, corresponding to the scaling regime of stratified turbulence identified in Billant 


& Ghomaz (2001) (see also (Lindborg 2006)). Another interesting approach consists of 


computing a full frequency-wavenumber energy spectrum in the four dimensional Fourier 
space (Glark di Leoni & Mininni (2015); see also|Gampagne et al. (2015) for an analogous 
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analysis in a laboratory experiment of rotating flow). This allows to directly assess how 
the energy is partitioned between modes which exhibit phase coherence, like waves, and 
other modes. However, this is a computationally expensive technique which requires a 
good time resolution in addition to a good space resolution (and hence, a lot of disk 
space). 

In this context, a pioneering analytical and numerical study of strongly rotating strat¬ 
ified turbulence in terms of normal modes was done by |BarteIIo (1995). The link between 
geostrophic adjustment at large scale and the inverse cascade of energy due to the conser¬ 
vation of (linearized) potential vorticity, is clearly established through a decoupling of the 
(fast) wave modes and vortical (slow) modes. By truncating the nonlinear interactions to 
some specific subsets of modes, as for example between three slow modes that lead to the 
inverse cascade of energy ( Charne^|1971 ), dominant effects can thus be identified. It is 
shown in particular that there are nonlinear interactions using a slow mode as a catalyst 
to push the wave energy to small scales (whereas triads involving three fast modes are 
negligible); furthermore, as the rotation is increased, the energetic exchanges between the 
slow vortical modes and the fast wave modes are weakened and such modes can be seen 


as actually being progressively decoupled. It was noticed by Smith & Waleffe (2002) that. 


for long times, the growth of large-scale energy is associated with geostrophic modes for 
1/2 < iV// < 2 (in which domain three-wave resonances vanish completely), whereas for 
large N/f, the ageostrophic thermal winds, also called vertically sheared horizontal flows 
(VSHF), are dominant, with no variation in the horizontal and no vertical velocity. One 
should note that if the forcing wavenumber is not well separated from the box overall di¬ 
mension, finite-size effects can develop that quench the possible (discrete) resonances, as 
already noted in the context of water-wave turbulence (see e.g. [Kartashova et al. (2008) 
and references therein). 


In the purely stratified case, it was found by Waite & Bartello (2004 20066) that 


the vertical spectrum of the waves moves to larger vertical wavenumber as stratification 


increases, in agreement with the argument found by Billant & Chomaz (2001) that the 
characteristic scale in the vertical is proportional to U/N, with U a typical large-scale 
velocity; at smaller scales the structures break-down in a nonlinear fashion although 
isotropy may not be restored, depending on the buoyancy Reynolds number and whether 


the Ozmidov scale is resolved (see § 2.2 below for definitions): the buoyancy scale Lb = 
U/N must be sufficiently resolved for over-turning to take place. The case of comparable 
(but not identical) strengths for rotation and stratification was analyzed by [Sukhatine 
I&: Smith (2008), where it is shown that there is a transition in behavior for N/f « 1: for 
smaller values, large scales are dominated by the vortical mode and small scales by wave 
modes, whereas for larger values of N/f, the wave modes dominate at all scales smaller 
than the forcing. 

The scaling behavior of slow and fast modes for rotating stratified flows is analyzed 
Kurien & Smith (2012[ 2014) for regimes where potential vorticity is linear, for small 


aspect ratio of the domain size in the vertical (6 = L||/L_l = 1/4), but with a unit 
Burger number Bu = SN/f, using high resolution direct numerical simulations and 
hyper-viscosity. When compared to similar flows with 6 = 1, differences arise, as for 
example when examining the steepness of the wave component of the energy spectrum, 
viewed as a signature of the layered structure of the flow. The case 6^1 was also 


examined for either purely rotating (Deusebio et al. 2014) or purely stratified flows (Sozza 
et al. ]M^ , where it is shown that the aspect ratio can affect the strength of the inverse 
cascade of energy to large scales. The combined effect of varying N/f and d on the 


formation of characteristic structures was further studied in Kurien & Smith (2014) for 
low Froude number (Fr = 0.002, ensuring that potential vorticity is again linear, i.e. 
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for low effective turbulence) and varying either the Rossby number at fixed 6 or vice- 
versa (in which case, the choice was made of N=f). The layering of the flow is clearly 
dependent on the strength of the rotation (Pr being fixed), and the wave and slow-mode 
spectra are quite different (when the total spectra are more difficult to distinguish). The 
inverse cascade of energy is clearly associated with the vortical (slow) mode, as shown 
as well in Brunner-Suzuki et al. (20141 for a variety of forcing functions in the context 
of sub-mesoscale oceanic mixing. It is also shown in Brunner-Suzuki et al. (2014) that 
small-scale (shear) layers do not directly modify the inverse cascade although it could 
strengthen it when bringing together nearby vortices through nonlinear advection when 
they are sufficiently numerous. 

Our study uses similar methods as the above mentioned authors, but pursues a slightly 
different goal: here, we focus on the possibility of an inverse cascade and its relation with 
the inviscid invariants of the system. We are interested in understanding the difference 
between the rotating-stratified case, for which an inverse cascade is possible ([Marino et al. 


2013), and the purely stratified case, where there is no such thing (although there is some 


transfer of energy towards the large scale). This difference has already been considered 
from the point of view of anisotropic energy transfers (Marino et al. ]M4t ; here we adopt 
a different point of view and we investigate the role of vortices and waves in the inverse 
cascade (or lack thereof). We also aim at characterizing the transition between these 
two regimes. Indeed, although previous studies have established that when an inverse 
cascade is present, it is due to the slow vortices, it is not clear how this phenomenology 
breaks up. We shall see that in the regime of parameters considered, we can identify 
three potential candidates, which can also be superimposed. To do so, we examine data 
stemming from direct numerical simulations of rotating stratified turbulence for a variety 
of N/ f ratios (and Froude and Rossby numbers) in cubic boxes at a Reynolds number 
of 1000. After describing the relevant equations and the numerical simulations in § [^ 
we review the normal mode decomposition in § [^ and analyze the numerical data in 
terms of the normal modes in §[^ In §[^ we discuss the results, first in the framework of 
the inviscid invariants of the system, and then by considering the effect of the buoyancy 
Reynolds number on the energy partition between slow and fast modes. Finally, our 
conclusions are exposed in § [^ 


2. Numerical Data 

2.1. Equations of motion and direct numerical simulations 

We consider incompressible flows in a reference frame rotating with angular velocity Q 
around axis (we note D = and / = 217 the Coriolis parameter). The gravity 
field is anti-aligned with the direction of rotation: g = —ge^. We work in the Boussinesq 
approximation: the density is assumed constant, equal to po > 0, except in the buoyancy 
force term, and we assume a linear background stratification: the density field is given 
by p(x) = Po N the Brunt-Vaisala frequency. Note that we are 

only considering the case of stably stratified flows. We introduce the temperature field 
(with the dimension of a velocity) 0(x) = The velocity field is denoted by u. 

Given the assumptions mentioned above, the equations governing the dynamics of the 
fields u,0 are as follows: 

9tu + u • Vu = —VP -I- z/Au — 2176^ x u — NOe^ + F, 
dtO + u • = Nuz + kAO, 

V • u = 0, 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 
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Id 

rip 

/ 

N 

Fr 

Ro 

Bu 

Re 

TZb 

Rui 

kp 

Rl/4 

512® 

44 

11 

0.32 

0.08 

0.25 

1000 

102.4 

2.53 

22 

Rl/2 

512® 

44 

22 

0.16 

0.08 

0.5 

1000 

25.6 

2.53 

22 

R3/2 

512® 

44 

66 

0.05 

0.08 

1.5 

1000 

2.5 

2.53 

22 

R3 

512® 

44 

132 

0.027 

0.08 

3 

1000 

0.73 

2.53 

22 

R4 

512® 

44 

178 

0.02 

0.08 

4 

1000 

0.4 

2.53 

22 

Fl/2 

512® 

178 

89 

0.04 

0.02 

0.5 

1000 

1.6 

0.63 

22 

F3/2 

512® 

58 

89 

0.04 

0.06 

1.5 

1000 

1.6 

1.90 

22 

F2 

512® 

44 

89 

0.04 

0.08 

2 

1000 

1.6 

2.53 

22 

F3 

512® 

29 

89 

0.04 

0.12 

3 

1000 

1.6 

3.80 

22 

F4 

512® 

22 

89 

0.04 

0.16 

4 

1000 

1.6 

5.06 

22 

F7 

512® 

13 

89 

0.04 

0.28 

7 

1000 

1.6 

8.85 

22 

FIO 

512® 

9 

89 

0.04 

0.4 

10 

1000 

1.6 

12.65 

22 

F20 

512® 

4.5 

89 

0.04 

0.8 

20 

1000 

1.6 

25.3 

22 

Foo 

512® 

0 

89 

0.04 

OO 

CX3 

1000 

1.6 

OO 

22 


Table 1. Model run parameters; names of runs at fixed Rossby number {resp. fixed Froude 
number) begin with a R {resp. F), followed by its N/f value. Up is the number of grid points, 
N and f are the Brunt-Vaisala and the inertial frequency (with / = 2fl, being the imposed 
rotation), Fr = U/[LN] and Ro = U/[Lf] are the Froude and Rossby numbers based on 
large-scale data, and the Burger number is defined as Bu = N/f for a cubic box as used here; 
TZb = ReFr^ and R^j — Re^^^Ro — ujrms/f are, respectively, the buoyancy Reynolds number 
and the micro-Rossby number, and kp is the forcing wavenumber. The runs belonging to the 
no-resonance zone have a shaded background. 


where P denotes the (rescaled) pressure field, fixed by the incompressibility condition, 
1 / the viscosity and k the thermal difFusivity; we work with Prandtl number equal to 
unity: Pr = v /k = 1. Finally, F is an isotropic forcing introduced only in the momentum 
equation; it has a fixed amplitude in a prescribed narrow shell centered on wavenumber 


kp, and random but constant in time phases. Note that, contrary to other studies (Kurien 
fc Smith|2012 2014, for instance), there is no forcing in the buoyancy equation, but the 


forcing has a vertical component (unlike, e.g. Waite & Bartello (2004) or Sozza et al. 


(2015)) and does not satisfy balance relations, even at the lowest order (unlike e.g. Waite 


& Bartello (20066) where the forcing is only in the slow modes). The projection of the 
forcing on the waves and slow modes (see §[^ has the same magnitude; hence we are not 
favoring any of these modes a priori. 

We integrate equations ( 2.1|2.3 ) numerically with a pseudo-spectral code, the Geo¬ 
physical High-Order Suite for Turbulence (GHOST) in a tri-periodic cubic domain of 
length L = 2tt. GHOST is parallelized using an hybrid MPI/OpenMP method and scales 
linearly up to 10® compute cores on grids of up to 6144® points (Mininni et al. 2011). 
No model is used for sub-grid scales, and dissipation occurs through standard Laplacian 
terms. The runs considered here all have a resolution of 512®, with forcing at wave num¬ 
ber kp S [22,23], and therefore Reynolds numbers Re « 10®. The Rossby and Froude 
numbers vary as indicated in Table the runs can be grouped in two series, one with 
constant Froude number {Fr = 0.04) and varying N/f (and therefore Ro) and one with 
constant Rossby number {Ro = 0.08) and varying N/f (and therefore Fr). The isotropic 
energy spectra and fluxes of the above runs have already been analyzed in [Marino et al. 


(2013), providing evidence for the existence of an inverse cascade of energy in the pres¬ 


ence of rotation and stratification, but not in the purely stratified case (see also Marino 


et al. (2014)). 
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2.2. Characteristic scales 

Rotating stratified turbulence covers a wide range of physical regimes, which can be 
characterized by three non-dimensional numbers, besides Pr. Given a characteristic 
length scale L and a characteristic velocity U at that scale, we introduce as usual the 
Reynolds number Re = ULjv, the Rossby number Ro = Uj (fL) and the Froude number 
Fr = U/{NL). The Rossby and Froude numbers quantify the strength of the Coriolis 
and buoyancy forces, respectively. They correspond to the ratio between a wave period 
(/“^ and N~^, respectively) and the eddy turnover time L/U. 

A difficulty in rotating-stratified flows arises from the variety of characteristic scales 
one can construct. To start with, the Froude number can be based on a characteristic 
horizontal length scale or a vertical length scale. It is known that stratihed turbulence 
develops strong gradients in the vertical (leading to the ubiquitous layered structure 


of the flow, see e.g. 

Herring & Metais 

(1989)), in such a way that the vertical Froude 

number remains of order unity ( 

Billant & Chomaz 2001 

Lindborg 

2006), even when 


the horizontal Froude number is small. The corresponding vertical length scale is the 
buoyancy scale. Lb = U/N and can be identified readily in such flows (see e.g. 


Rorai 


et al. (2014)). In the presence of rotation, the layers are expected to be slanted (with 
respect to the horizontal). In that case, it is conjectured by Billant & Chomaz (2001) 
that the characteristic length scale in the vertical also involves both N and /, like a 
deformation radius (for which the effects of rotation and stratihcation balance each other, 
with Nk± ~ /fc||)- It was shown by Rosenberg et al. (2015) that a Froude number based 
on a vertical Taylor scale involving vertical velocity gradients, iz = \/ (u^)/{{dzU±_)‘^) 
(where the angular brackets denote spatial averaging) is also very close to unity. Here, 
we dehne the Froude and Rossby numbers (as well as the Reynolds number) based on 
the scale of the forcing, Lp = ‘1 'k jkp. Since the forcing is isotropic, there is no reason to 
distinguish between vertical and horizontal Froude number with this definition, although 
anisotropies and different characteristic scales in the horizontal and the vertical will of 
course develop spontaneously in the flow. It follows that for the purpose of our study, 
at unit aspect ratio, the ratio of the Rossby and Froude numbers Bu = Ro/Fr, referred 
to as the Burger number, reduces to the ratio of the Brunt-Vaisala frequency and the 
Coriolis parameter: Bu = N/f. 

The recovery of isotropy should occur at a scale for which the eddy turn-over time 
and the wave period become comparable. This allows one to define an Ozmidov and a 
Zeman scale ( |Zeman|1994 Mininni et a/.|20l2 ) written, respectively, for purely stratified 
or purely rotating flows, assuming an isotropic Kolmogorv range, Ev{k) 
small scale, as: 


at 



with cv = the kinetic energy dissipation rate. The Froude and Rossby numbers 

measure the ratio of these scales to the characteristic scale: Fr = (iozjL')^P,Ro = 
{izejL'i^P. Finally, if one wants to construct a measure of the actual development of 
small-scales in the flow, a buoyancy Reynolds number Hb and a micro Rossby number 
R,^ can be defined in the following way: 

Hb = ReFr^, R^ = = Re^PRo, (2.5) 

so that R^ is the equivalent, for rotation, of the buoyancy Reynolds number. 

With these definitions, and recalling the Kolmogorov dissipation scale = {P jcvY^^, 
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J_^^_L 

0 10 20 30 


t 
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-^N/f = 3 -^N/f = 4 


Figure 1. Time evolution of the potential enstrophy for the series of runs at constant Froude 
number (left) and at constant Rossby number (right) with various N/f. The vertical scale is in 
logarithmic coordinates for the constant Rossby number series. 


one can show that TZb = ■ In other words, when one wants to measure the 

intensity of stratified turbulence, the Ozmidov scale plays the role of the integral scale and 
TZb that of the Reynolds number. In particular, when TZb is of order one, like in most of 
the runs shown here (see Table[^, the Ozmidov scale coincides with the dissipation scale. 
Also note that, when « 1, the small-scale turbulence is at least as vigorous as the 
imposed rotation frequency (see, e.g., Sagaut & Gambon (2008)). These non-dimensional 
numbers are computed in Table 

One of the goals of this paper is to study how the respective role of waves and vortices 
(see below) depends on the ratio N/f. Note that there is a range that can be expected to 
have a special behavior: when | < y < 2, there are no three-wave resonances (the simple 
proof relies on the fact that the inertia-gravity wave frequency cr(k), defined by 3.4, is 
bounded by / and N] see for instance Smith fc Waleffe| (2002), § 6.1). Hence the no¬ 
resonance zone separates two regimes: N/f <1/2 and N/f > 2. Here we shall focus on 
the latter range. 


2.3. Potential Vorticity 
We define potential vorticity as usual: 

n = fdz9 — Nujz + oj ■ VO, 


( 2 . 6 ) 


where a; = V x u is the standard vorticity. Potential vorticity is a Lagrangian invariant 
of ideal Boussinesq flows: in the absence of forcing and dissipation, it is conserved along 
the trajectory of a fluid parcel, as expressed by the equation dWT -I- u • VH = 0. Hence, 
it is analogous to vertical vorticity in 2D turbulence, a fundamental difference being 
that potential vorticity is not a linear functional of the dynamical fields, because of the 
quadratic term II2 = u> ■ VO, in addition to the linear term Hi = fdzO — Nujz- It follows 
that, unlike the 2D case, the L 2 norm of potential vorticity, referred to as potential 
enstrophy and denoted by T, is quartic and is not necessarily conserved by the Fourier 
truncation (see, e.g. Aluie & Kurien (2011)). We show in Fig. the time evolution of 
potential enstrophy in the two series of runs presented here. In the series of runs at 
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-^N/f = 2 -^N/f = S -^N/f = A 
-^N/f = 7 -^N/f = 10 N/f = 20 
- N/f = oo 


Rosenberg and A. Pouquet 


Ro = 0.08 



t 


- N/f =1/A—-N/f = 1/2---- N/f = 3/2 

-^N/f = 3 -^N/f = 4 


Figure 2. Time evolution of the ratio of quadratic potential enstrophy to total potential en- 
strophy, for the series of runs at constant Fronde number (left) and at constant Rossby number 
(right) with various N/f. 


constant Fronde number, potential enstrophy decays to about 25% of its initial value, 
then seems to remain constant in time. The decay is faster for higher N/f ratios, but 
the final potential enstrophy level is about the same for all the runs. On the contrary, 
the series of runs at constant Rossby number exhibits very different levels of potential 
enstrophy, across about three orders of magnitude, and it increases with increasing N/f, 
i.e. decreasing Froude number. We conclude from these two graphs that, in the inverse 
cascade regime studied here, the value of potential enstrophy at late times is a function of 
the Froude number only. Also, potential enstrophy conservation is better satisfied when 
stratification is strong. We now introduce the quadratic, cubic and quartic components 
of potential enstrophy, dehned respectively as 



(2.7) 

r3 = 2 fniU,. 

( 2 . 8 ) 

r, = / nl. 

(2.9) 


such that F = r 2 + F 3 + F 4 . r 2 . Fa and r 4 are respectively quadratic, cubic and quartic 
in the field variables, but they all have the same physical dimension, that of frequency 
square or square vorticity. The ratio r 2 /r is shown in Fig. Note that this ratio can 
exceed unity, since the cubic component of total potential enstrophy is not positive 
definite (it measures correlations between the linear part of potential vorticity and the 
nonlinear part). This hgure shows that in all the runs, total potential enstrophy is well 
approximated by its quadratic part, in agreement with the findings of |Kurien &: Smith] 


(2014), although it is less and less the case as rotation decreases (constant Froude series of 
runs) or in the N/f = 1/4 (weak stratification) case, which also has the highest buoyancy 
Reynolds number « 100 (see Table [^. 

In 2D turbulence, enstrophy conservation has very important consequences: it is ul¬ 
timately responsible for the inverse cascade. It is a legitimate question to ask whether 
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potential enstrophy conservation places such a strong constraint on rotating-stratified 
flows. Another important aim of this paper is to show that this constraint is strongly 
dependent on whether there is rotation or not. 


3. Waves and Vortices: Normal Mode Decomposition 

3.1. The linear dynamics 

Inertia-gravity wave dynamics has been studied for a long time (see for example Hassel- 
(19661; Garrett & Munk (1979); Muller et al. (1986); Bartello (1995); Staquet & 


mann 


Sommeria (2002)) and is summarized very briefly here in order to recall essential concepts 
and to define notations. Indeed, the linear terms introduced by the Coriolis force and the 
buoyancy force in equations ( 2. Ip. 3 ) allow for the propagation of dispersive waves. The 
linearized inviscid and unforced dynamics reads: 


dtu = —VP — fez XU — NObz 
dtO = Nuz, 

V • u = 0. 


(3.1) 

(3.2) 

(3.3) 


This set of equations corresponds to the low Froude number limit (Pr —> 0) of the 
Boussinesq equations, if we assume that the time scale is the buoyancy time scale 
N~^ (Lilly 1983). The above equations admit travelling wave solutions of the form 
Uj(x,t) = = 0 ehkx-a(k)t)^p ^ p^g*(k.x-<T(k)t)^ for arbitrary wave 

vector k, with relations between the complex amplitudes Uj,Q,Po fixed by the above 
equations, and with the dispersion relation 


cr(k) = ^ pkl/k^ + N^kl/k^. 


(3.4) 


Such waves are called inertia-gravity waves (and include, of course, waves traveling in the 
—k direction). Note that although the Coriolis parameter only appears in the evolution 
equation for the horizontal velocity, and the Brunt-Vaisala frequency only for the vertical 
velocity and temperature, inertia-gravity waves couple the four dynamical fields because 
of the incompressibility condition. Without pressure, we would simply have inertial os¬ 
cillations, at frequency /, for the horizontal velocity, decoupled from gravity waves of 
frequency N acting on the vertical velocity and the temperature. 

The nonlinear term introduces a coupling between waves with different wave vectors. 
A possible approach is to treat the dynamical fields as a wave field with slow nonlinear 
evolution of the magnitude of each wave mode, using a perturbative approach. This is a 
first kind of turbulence, referred to as wave turbulence ( Zakharov et aL|1992 Nazarenko 
|2010 Newell fc Rumpf|[2011 ). 

Now, the set of linearized equations also admits stationary solutions, where the right 
hand side of equations ( 3.1|3.2 ) vanishes. These solutions have no vertical velocity {uz = 
0), and the pressure gradient balances at the same time the Coriolis force {geostrophic 
balance, V_lP = —fez x u) and the buoyancy force [hydrostatic balance, dzP = —N6ez). 
In other words, rotation and stratification act to maintain horizontal motion. This kind 
of balanced motion can also be seen as a low Froude number limit [Fr —> 0) solution 
of the Boussinesq equations, assuming that Bu = 1 and that the time scale is the eddy 


turnover time (Lilly 1983). It corresponds to eddies of typical aspect ratio f/N, which 


are advected by the horizontal velocity field, much like vortices in 2D turbulence. Hence, 
another kind of turbulence, referred to as geostrophic turbulence, is encapsulated in the 
Boussinesq equations. 
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3.2. The normal mode decomposition in Fourier space 

The two kinds of motion, waves and vortices, described above, correspond to normal 
modes of the linearized Boussinesq equations. Although they are coupled in the nonlinear 
dynamics, they provide a useful framework for analyzing the full system. Here, we shall 
try to disentangle the behavior of these two contributions in the simulated fields. To 


separate the two contributions, it is convenient to express them in Fourier space (Bartello 


1995 


see also Herbert et al. (2014)): we introduce B = lA' the set of wave vectors, and 


for each wave vector B, the vectors 


Zo(k) =M(k) 


Z_(k) = M(k) 


and Z+(k) =M(k) 


(3.5) 


where the expression of the matrix M(k) is given by: 


M(k) = < 


/ Jk 2 k\\-ikik\\<T{\i.) 
-/feifcii - ifc2fc||cr(k) 


\/ 2 kk± 'T(k) 


1 


1 

0 

Vo 


V 

0 

0 

0 


zfc^cr(k) 

-Nk'i 


-V2Nk2k± 

'/2Nkik± 

0 

-V^fk\\kj_ 


/fc2fc|| + ifcifc||tT(k) ^ 
—/fclfc|| + lfc2fc||0'(k) 

—zfc^cr(k) 


-Nkl 




1 

0 

Oj 


if fc_L ^ 0 


if fc_L = 0 


(3.6) 


The vectors Zo(k), Z_(k) and Z-|_(k) are mutually orthogonal, and normalized to unity: 

Z,(k)tZ,(k) =(5„, (3.7) 

where f denotes the transpose and complex conjugation. Although the matrix M(k) is 
rectangular, it satisfies the hermitian identity M(k)^M(k) = I 3 . They are the normal 
modes of the linearized dynamics: their evolution in the linearized dynamics is given by 


L(k)Zo(k) = 0, 

L(k)Z+(k) = za(k)Z+(k), 
L(k)Z_(k) = -icr(k)Z_(k), 


(3.8) 

(3.9) 
(3.10) 


where L(k) is the linear operator in Fourier space associated to the linearized dynamics 
give by Eqs. 3.1|3.3 We see that the mode Zo(k) is conserved by the linear dynam¬ 
ics, while the modes Z±(k), corresponding to the inertia-gravity waves, propagate with 
frequency cr(k). Therefore, Zg, whose characteristic evolution time is the eddy turnover 
time, is referred to as the slow mode. The normal modes can be used as a basis for 
the Fourier modes of the dynamical fields. Introducing the Fourier decomposition of the 
velocity and temperature fields 


Mi(x) = ^ Mi(k)e*‘"'’‘, 
keB 

0(x) = ^ 0(k)e 
keB 


(3.11) 


^zk-: 


(3.12) 
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and denoting X(k) = *('Ui(k), U 2 (k), M 3 (k), 0(k)) the Fourier coefficients for wave vector 
k, they can be expressed as a linear combination of the normal modes: 

X(k) = Ao(k)Zo(k) + ^_(k)Z_(k) + A+(k)Z+(k), (3.13) 

M-(k)\ 

= M(k)Uo(k) . (3.14) 

\A+{k)J 

A(k) 

Introducing the projectors Po(k) = Zo(k)Zo(k)I andPw(k) = Z_(k)Z_(k)I+Z+(k)Z+(k)I, 
we have of course Po(k)X(k) = Ao(k)Zo(k) and Pw(k)X(k) = (k)Z_ (k)+A+(k)Z+(k), 

and the two projectors are mutually orthogonal, so that Po(k) 0 Pw(k) = I a- Denoting 
^ : {ui,U2 ,U3,9) I—>■ (ui,U2,U3,0), with x : k G B ^(k) G C, the Fourier transform, 
and #0 : {ui,U 2 , 0) (k Po(k)X(k)), : (^i, ^ 2 , U 3 , 6) ^ [k^ P, 4 /(k)X(k)), 

the operators o 3 ^^ o ^ and o o ^ are mutually orthogonal 

projectors acting on the set of dynamical fields: 

■^0 “ "^0 ) ~ (3.15) 

= 0 , 3^^®^w=ld- (3.16) 

Hence, we may write u = Uo0uw and Q = 9o-\-0w, with (ug, Ug, Uq, 9o) = 3^o(u^,u^, 9), 

and (ulf^,u^,u^,9w) = 3^w(u^,u^,u^,9). It is easily checked that the projection on 
the slow modes defined above has vanishing vertical velocity: Mg = 0, and is divergence 
free: V • Ug = Vj_ • Ug = 0. It follows from the above that the total energy 

Et = \ [{u^+9^) = ^J2 ^(k)^X(k) = i 5] A(k)t A(k) (3.17) 

d ^ keB ^ ksB 

breaks up in two pieces: Et = Eq -\- E\y, where Eq is the energy in the slow modes: 

Eo = l [iul + 9l) = ^ 5]|Ag(k)|^ = 1 ^ F;g(k), (3.18) 

d kGB ^ kGB 

and Ey^ the energy in the wave modes: 

Ew = l f (ul. + 9\) = i ^[|A_(k)|^ 0 |AH.(k)|^] = 2 E ■ (3-19) 

d ^ kGB ^ kGB 

Note that due to the orthogonality property of the projections, the cross terms in physical 
space cancel out: f Ug • uw + f 9o9iv = 0. However, these terms do not vanish individ¬ 
ually. In other words, the vector (ug,0g) is indeed orthogonal to the vector (Uu,,^^,), 
but this does not imply that the vector Ug is orthogonal to u^, where orthogonality 
refers to the canonical scalar product in each of these two Hilbert spaces. A consequence 
is that care should be taken when trying to mix decomposition of the energy in terms 
of kinetic/potential energy and in terms of vortical/wave modes. Here, we shall avoid 
altogether distinguishing between kinetic and potential energy. 


3.3. Normal mode decomposition, balance and potential vorticity 
As can be checked explicitly ( Herbert et fl/.|20'l4 e.g.), in the linear framework, the slow 
modes satisfy balance relations: in the presence of rotation, they are in both geostrophic 
and hydrostatic balance. By contrast, for purely stratihed flows, the slow modes are not 
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in hydrostatic balance, except for purely vertical wave vectors, i.e. for horizontally ho¬ 


mogeneous modes (the VSHF modes, see § 4.3). Of course, nonlinear effects break these 


relations, and a flow initialized in the image of the projector (i.e. containing no wave 
modes initially) will eventually lose balance and radiate inertia-gravity waves. If one is 
interested in maintaining the balance relations, as much as possible, without generating 
such oscillations, the initial state should contain a wave component designed to compen¬ 
sate the nonlinear effects; such nonlinear initialization procedures have been designed in 


numerical weather prediction (Baer & Tribbia 1977 Leith 1980). Alternatively, balanced 
models, where the dynamics of the waves is slaved to that of the slow modes, can be 
written to remain on a slow manifold in phase space ( ]Warn et a7||199'5 Vannest^|2013 ). 
For the sake of simplicity, we shall be content here with projecting the data on the vector 
space spanned by the slow modes (which can be seen geometrically as a tangent subspace 
for the slow manifold, or as a zeroth order approximation to the slow manifold) and its 
orthogonal subspace, spanned by the wave modes, as defined above. This means that 
although the normal modes introduced above are always well defined, care should be 
taken when interpreting the projections in terms of balance relations. 

The normal modes are also related to potential vorticity. An argument due in particular 
Smith & Waleffe (2002) relies on the fact that in the linear dynamics. Hi is a point- 


to 


wise invariant: djli = 0. It follows by applying the Fourier transform in both space and 
time, that (Tlli(k, cr) = 0, where cr is the time frequency. Hence, wave modes, for which 
(T(k) 7 ^ 0 must have Hi (k, a) = 0, while vortical modes (i.e., by definition, those for which 
ni(k, cr) 7 ^ 0) must have zero frequency. In fact, it is easi ly showri that in general, the 


identity IIi(k, f) = —zA:cr(k)Ao(k, t) holds, and therefore (Bartello 

Mil): 


1995 


Herbert et al. 


= fc^cr(k)^|Ao(k)|^ 


(3.20) 


keB 


It follows that, even when considering the full (nonlinear) dynamics, waves do not carry 
any (linear) potential vorticity. Similarly, we know that in the presence of rotation, all 
the slow modes contribute to (linear) potential vorticity. In the purely stratified case, 
however, the slow component of the VSHF modes {k± = 0) contributes neither to (linear) 
potential vorticity, nor, therefore, to (quadratic) potential enstrophy r 2 . All the other 
slow modes still carry (linear) potential vorticity. 

The fact that the connection between slow modes and vortical modes (i.e. modes which 
contribute to linear potential vorticity) does not depend on whether we are considering 
the linear or the full nonlinear dynamics, while the connection between slow modes 
and balanced modes does, is due to the kinematic nature of the connection between 
linear potential vorticity and the dynamical fields. By contrast, the notion of balance is 
dynamical in the sense that it involves pressure, which is determined by the equations 
of motion. When including the nonlinear term in the dynamics, another component of 
pressure appears, which breaks the balance relations. 


4. Partition of energy between slow and wave modes 
4.1. Time evolution of global quantities 

The time evolution of the two energy components, wave and vortical, is shown in Fig. 
(see also Fig. for the series of runs at constant Froude number (left) and at con¬ 
stant Rossby number (right), excluding the runs for which the three-wave resonances 
are suppressed. Although the initial conditions all have about the same energy in each 
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Fr = 0.04 



- Eo - Eyy N/ f = 2 

-^N/f = 3 ^N/f = 4 -^N/f = 7 
-^N/f = 10 N/f = 20 - N/f = oo 


Ro = 0.08 



t 


- Eo ---Ew - N/f = 1/4 

-~N/f = 3^N/f = 4 


Figure 3. Time evolution of the wave (dashed lines) and vortical (solid lines) parts of the total 
energy for the series of runs at constant Fronde number (left) and at constant Rossby number 
(right) with various N/f. Runs within the no-resonance zone have been excluded for clarity. 
Note that the curves corresponding to the decay of the wave energy almost collapse onto a 
single curve: see Fig. for a figure in logarithmic coordinates where they can be distinguished. 


■ 10 “^ 




Figure 4. Left: Time-averaged growth rate of the energy in the slow modes as a function of 
N/f, for the two series of runs, at constant Froude number (circles) and at constant Rossby 
number (squares), respectively. Right: Time averaged fraction of energy in the slow modes as a 
function of N/f, for the two series of runs, at constant Froude number (circles) and at constant 
Rossby number (squares), respectively. The no-resonance zone is the shaded area. 


component, we see that in all the runs, the wave energy decays very fast and the slow 
modes dominate, even for the runs that do not undergo a strong inverse cascade (in the 
purely stratified case, there is no inverse cascade at all). This is also true for the runs in 
the series at constant Rossby number (Fig. right) and in the no-resonance zone (not 
shown). After an initial phase, we observe in the runs for which N/f is small enough 
{N/f < 10) a linear growth of the energy in the slow modes, characteristic of an inverse 
cascade. For N/f > 2, the growth rate appears to be a decreasing function of N/f and 
this is confirmed by plotting the time average of the time derivative of the energy in 
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Fr = 0.04 Ro = 0.08 
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Figure 5. Time evolution of the wave component of the total energy in logarithmic coordinates, 
for the series of runs at constant Froude number (left) and at constant Rossby number (right) 
with various N/ f. 


the slow modes (Fig. left). This figure also confirms that the growth rate vanishes at 
N/f = 20 or in the purely stratified case (see also Marino et al. (2014)). There is also 
a linear growth of the energy in the slow modes for the runs in the no-resonance zone 
in the series at constant Froude and for all the runs in the series at constant Rossby. 
For both series of runs, the growth rate appears to reach a maximum located inside the 
no-resonance zone (Fig.left). In addition, the growth rate seems to depend primarily 
on N/f rather than the Froude and Rossby numbers. These conclusions are consistent 


with the analysis of the time evolution of the kinetic energy obtained earlier (Marino 
et qZ.|[ 20T^. Figure]^ (right) summarizes how the fraction of energy in the slow modes. 


averaged in time (over all the snapshots following the peak of enstrophy), evolves as a 
function of N/ f. Surprisingly, the case where this fraction attains its minimum value, 
although it is still about 0.9, is not the purely stratified case, but the rotation domi¬ 
nated case N/f = 1/4, which also has the highest TZb)- Since Fig. [^does not allow to 
distinguish the different curves for the wave energy at various N/ f ratios, we show this 
quantity in logarithmic scale in Fig.[^ For values of N/ f above 7 (this is also true for the 
case N // = 1 /4) , we see a decay phase followed by a regime where the level of energy in 
the waves stays constant. For the other N/f values, the wave energy drops suddenly at 
the end of the decay phase, then increases back to more or less its value before the drop, 
and seems to stay constant afterwards. Although the time resolution is not sufficient for 
precise measurements, the decay seems to follow roughly a power-law, and it is steeper 
for smaller N/ f ratios (except for the N/ f = 1/4 and N/f = 1/2 cases, at substantially 
higher TZb)- 


4.2. Spectral analysis 

We show in Fig. [^the isotropic spectra of energy in the wave and slow modes, denoted 
respectively E\\r{k) and E^Qt), for the series of runs at constant Froude number. The 
isotropic spectra are defined, as usual, as integrals over spherical shells in Fourier space: 

Eo{k)= f |Ao(k)|2dk, Ewik)= f [|A+(k)|V|A_(k)|"]dk. (4.1) 

Jj|k||=fc "'l|ki|=fe 
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Figure 6. Isotropic energy spectra of the wave (solid lines) and slow mode (dashed lines) for 
the series of runs with constant Fronde number {Fr = 0.04) and varying N/ f. The times shown 
are indicated in the legend at the bottom, going from purple to blue, green, yellow and red in 
increasing order. Two different spectral laws are indicated as well for comparison. 


For all the runs, at early times {t = 1.5, before the peak of enstrophy), the slow mode 
component dominates at large scales while the wave component dominates at small scales. 
Note that the wavenumber at which the two spectra cross increases with N/f. When the 
turbulence develops, the wave energy spectrum decays rapidly, including at small scales, 
and the slow mode spectrum dominates at all scales. However, the wave energy at small 
scales increases when N/f increases. The other salient feature of the isotropic energy 
spectra is the growth of the slow mode energy at scales larger than the forcing scale. 
In most of the runs, up to about N/f = 7, this growth takes the form of an inverse 
cascade, with a spectral slope that is consistent with However, the inverse cascade 

of the slow modes weakens as N/ f increases, as already found from the analysis of global 
quantities. Already, at N/ f = 7, the spectrum at the final time of the simulation has a 
narrow range, and at N/f = 10, the energy in the slow modes seems to saturate 

at large scale with a flat spectrum k^. Note that we cannot rule out completely the 
possibility that the —5/3 range would extend further towards large scales with a larger 
integration time. However, it should be noted that the growth rate of the energy in the 
slow modes for these values of N/f is very small (see Fig. and that the slow mode 
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energy spectrum have not evolved much between t = 22.5 and t = 30. This indicates 
that even if the —5/3 range was to extend towards larger scales, it would require a 
very long integration time. Besides, the growth does not seem to be self-similar, since 
the spectral slope at scales larger than the —5/3 range becomes shallower (almost flat), 
unlike traditional inverse cascade phenomena. The flat spectrum for slow mode energy 
at large scales is the only notable feature of the N/ f = 20 run, and it is also present in 
the purely stratified run, although in that case there is an accumulation of energy in the 
k = 1 modes. Such a flat spectrum was already observed for kinetic energy in Smith 


20141). 


Waleffe ( |2002 ) for N/f = 10 (see also Kimura & Herring (2012) and Marino et al. (2013 


The transition between these two scenarii is further described in Fig. which shows 
the spectral slope a of the isotropic energy spectrum of the slow modes at large scale, 
defined by Eglk) ^ fc“ for k < kp, and computed numerically by a least-square fit, as a 
function of N/ f. This figure shows that the transition appears to be continuous, although 
the specific values of the spectral slopes lack in accuracy since the inertial range is rather 
short (remember that kp G [22,23]), and since they may depend on the details of the 
least-square fitting procedure: choice of snapshot, excluding the fc = 1 modes to get rid 
of the formation of a possible condensate at the largest scale,... In particular, no time 
averaging has been performed here. Higher resolutions would also provide improvements 
in accuracy. Another caveat is that, as noted above, the intermediate values (e.g. JV/f = 7 
or 10) may move closer to —5/3 with a larger integration time. 

The inverse cascade of the slow mode energy spectrum is also seen in the series of runs 
at constant Rossby number (Fig. [^. Interestingly, for the runs with small N/ f ratios 
{N/f = 1/4 or 1/2), i.e. the runs with the highest Froude numbers, the waves dominate 
at small scales even at late times, when these scales seem to have attained stationarity. 
This is probably a consequence of the higher buoyancy Reynolds number, rather than 


the Burger number (see the discussion in § 5.3). 


In addition to the isotropic spectrum, obtained by summing over spherical shells in 
Fourier space, it is always relevant for anisotropic flows to compute the so-called reduced 


urier space (see e.g. 

Clodeferd & Staquet 

2003 

Sagaut & Cambon 

2008 Mininni 

Mininni et a/.||2012 

Marino et al. 

2014 

. Introducing the axisymmetric spectra 


eo{kj_,k\\),ewik±,k\\) defined by 


(4.2) 


r 2 

eo{kj_,k\\) = / |Ao(A:_l cos^e^; + A:_l sin(/)ey + A:||e_3)| k±d(j)^ 

Jo 

p27r 

ewik±,k\\)= / [\A^{k± cos (pSx k±sin(l)ey k\\ez)\ \A-{k± cos (j)exk±sm(j)eyk\\ez)\ ]k±d(j), 

Jo 

(4.3) 


the perpendicular and parallel spectra are defined as follows for each type of mode, 
s = 0,W: 

p^max 

Es{k±) = / es(fc_L, fc||)dfc||, Es{ki\) = / es{k±,ki\)dk±. (4.4) 

Jo Jo 

Figure shows the perpendicular spectrum of energy in the wave {Ewikp)) and slow 
modes {EQ{k±)); again, except for the runs with the highest Froude number {Fr = 0.32 
and Fr = 0.16, not shown), and therefore, buoyancy Reynolds number, the slow modes 
dominate at all horizontal scales after the peak of enstrophy. In all cases, the perpen¬ 
dicular spectrum of the slow modes undergoes a growth of energy in horizontal scales 
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N/f = lli,Fr = f).‘i2 JV//= 1/2, Fr = 0.16 



Figure 7. Isotropic spectra of the wave and slow mode components of energy for the series of 
runs with constant Rossby number [Ro = 0.08) and varying N/f. The symbols are the same as 
in Fig. 1^ 


larger than the forcing scale (note, however, that since the forcing acts on a spherical 
shell in Fourier space, there is some injection of energy at all the horizontal scales such 
that 0 < k_\_ < kp). As for the isotropic spectrum, this growth appears to be self-similar 
for N/ f ratios smaller than 10. As N/f increases, the perpendicular spectrum saturates 
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AT// = 1/4, ijo = 0.08, Fr = 0.32 



N/f = 7 ,Ro = 0.28, Fr = 0.04 



Rosenberg and A. Pouquet 


N/f = 2,Ro = 0.08, Fr = 0.04 



N/f = oo,Ro = 00 , Fr = 0.04 



- Eo--- Ew -1 = 1.5- 1 = 7.5 - 1= 15-1 = 22.5-i = 30 


Figure 8. Perpendicular spectra of the wave and slow mode components of energy for selected 
values of the parameters. Scaling laws are represented only for reference. Line encoding and 
colors are similar to the two preceding figures. 


at large horizontal scales; see the purely stratified case in Fig.j^for instance. The spectral 
slope in this range varies depending on the non-dimensional parameters, and it is difficult 
to provide precise values here. Nevertheless, it appears that the perpendicular spectrum 
becomes shallower as N/f increases. The spectral index aj_, defined by Eo{k±) ^ 
for k± < kp, and computed numerically by a least-square fit, is represented as a function 
of N/f in Fig. 10 The caveats mentioned above regarding accuracy of this spectral index 
still apply here. For instance a slightly steeper spectral slope might have been obtained in 
the purely stratified case if we had chosen a range of wave numbers closer to the forcing 
wavenumber: Marino et al.\ ( |2014 ) claim that the kinetic energy (here we are showing 
the slow mode energy) has a —5/3 slope in this range. The perpendicular wave spectrum 
differs markedly in the rotation-dominated, weakly stratified case and in the strongly 
stratified cases. In the former case, more energy remains at small horizontal scales, while 
in the latter case it falls off sharply at scales smaller than the forcing scales, and it is 
shallow or flat at scales larger than the forcing scales. 

Similarly, the parallel spectra Ewik\\), Eo{k\\) of energy in the wave and slow modes, 
respectively, are shown in Fig. Again, once the turbulence sets in, the slow modes 











Waves, vortices, and the inverse cascade of rotating-stratified turbulence 19 
N/f = 1/4, _Ro = O.OS, Fr = 0.32 N/f = 2.Ro = 0.08, Fr = 0.04 



N/f = 4, Ro = 0.16, Fr = 0.04 N/f = oo, Ro = oo,Fr = 0.04 



- Eo--- Ew -i = 1.5-1 = 7.5-1 = 15- 1 = 22.5-1 = 30 


Figure 9. Parallel spectra of the wave and slow mode components of energy for selected values 
of the parameters. Scaling laws are represented only for reference. Line encoding and colors are 
similar to the three preceding figures. 


dominate at all vertical scales except in the two cases with high buoyancy Reynolds 
numbers, where the waves dominate at small-scales (for the case N/f = 1/4, the waves 
dominate roughly at all the vertical scales smaller than the forcing scale kp, although 
the energy in these scales decay extremely fast). Overall, the energy in the small scales 
increases as we move from the rotation dominated case (where the spectrum is very 
steep) to the purely stratified case (note the different vertical scales on the panels of 
Fig.§. This statement also holds for the wave component of the energy: while the waves 
have most of their energy at large-scales for rotation dominated cases or cases where 
rotation and stratification are comparable (i.e. relatively low N/f ratios), they become 
mostly small-scale when stratification prevails. The buildup of slow mode energy at large 
vertical scales occurs differently from the perpendicular spectra: already at very early 
times, the parallel spectrum of the slow modes is more or less flat at scales larger than 
the forcing scale. Besides, for the lowest N/f ratio (1/4), there is relatively little change 
in the slow mode parallel spectrum at large scales. As for the perpendicular spectra, there 
does not seem to be a universal scaling law in this range: the spectral slope ay (defined 
by Eo{k\\) ^ for /cy < kp) varies between about —2 and —5/3 until N/f = 4 (see 
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—*— Q, Fr = 0.04 
—©— Q, Ro = 0.08 

- - Q;j_,i^r = 0.04 
_£>. q;_l, i?,0 = 0.08 

- - Q|l, Fr = 0.04 

- •4' • ■ Q;||, Ro = 0.08 


Figure 10. Spectral indices for the isotropic (solid lines), perpendicular (dashed lines) and 
parallel (dash-dotted lines) vortical energy spectra at scales larger than the forcing, as a function 
of N/f, for the two series of runs: constant Fronde number (star, x and -I- symbols), and constant 
Rossby number (circle, square and diamond symbols). The no-resonance zone is the shaded area. 



■ Fr = 0.04 
- Ro = 0.08 


Figure 11. Fraction of energy in the vertically sheared horizontal modes {k± = 0 modes) as a 
function of N/f, for the two series of runs. The no-resonance zone is the shaded area. 


Fig. 10), then the spectrum becomes much shallower at large scales, and it is almost flat 


for the purely stratified case and for the case with N/f = 20. 

4.3. Role of the VSHF modes 

Several studies of stratihed flows have underlined the role of horizontally homogenous 


modes: 

Herring & Metais (1989 

1 and 

Staquet & Godeferd 

(1998 

1 have reported a tendency 

towards layering, while Smith & Waleffe (2002 see also 

Godeferd & Gambon 

(1994 

) for 


an EDQNM perspective) have shown that at very long time, energy accumulates in the 
horizontally homogeneous modes (defined by k± =0), with a peak at some value of fey 
which defines the number of horizontal layers in the vertical, hence the name vertically 
sheared horizontal flows (VSHF). Our numerical simulations are stopped earlier in order 
to avoid the accumulation of energy in the gravest mode of the system that is inherent to 


all inverse cascade phenomena (see e.g. Boffetta & Ecke (20f2) for a recent review in 2D); 


this provides an incentive for looking at these modes in the direct numerical simulations 
carried out here. 

By definition, the VSHF modes are horizontally homogeneous, but may depend on the 
vertical coordinate. In general, they have an inertial wave component (whose frequency is 
the Coriolis frequency /) corresponding to an oscillation of the horizontal components of 
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velocity, with vanishing vertical velocity, and a slow mode component that is proportional 
to temperature fluctuations. Hence, their wave component is purely kinetic and their slow 
mode component is purely potential. When rotation vanishes, the dispersion relation also 
vanishes and there is no oscillatory component anymore. 


The ratio of energy in these VSHF modes is shown in Fig. 11 as a function of N/ f for 


the two series of runs. It is striking that it remains of the same order of magnitude (and 
relatively small; it is in fact negligible compared to the waves for instance, see Fig. 
for all finite values of fV//; it is only when we move to the strictly non-rotating case 
that the energy in the VSHF modes jumps by about two orders of magnitude. Then, 
the energy in the k± = 0 modes becomes comparable to the energy in waves with any 
other wave vector. Note that, although we do not have enough runs at high buoyancy 
Reynolds number for a systematic study, it does not seem to be a relevant parameter as 
far as this issue is concerned since the two runs with higher buoyancy Reynolds have a 
similar level of energy in the VSHF modes as the other runs. 

4.4. Energy fluxes 

Let us now turn to the fluxes of energy in the different modes. We first define the transfer 
functions ro(k) and Tw(k): 


To(k) =3? 


x„(k)t jPuj(p)X(q) 

p+q=k 


(4.5) 


ru.(k)=5ft Xu^(k)t ^ ik^uflp)Xiq) , (4.6) 

p+q=k 

where we have introduced the projections Xo(k) = Po(k)X(k) and X^^A(k) = Pw(k)X(k). 
The evolution of the energy in the slow and wave modes at wave vector k is simply gov¬ 
erned by 

{dt + 2,yk^)Eo{k) = To{k), (4.7) 

{dt + 2iyk^)Ew{k) = Tw{k). (4.8) 

As usual, we can define the isotropic, perpendicular and parallel transfer functions Ts(k), 
Ts{k±) and T's(fc||) (with s = 0,W) by integrating over a spherical shell, a plane or a 
cylinder in Fourier space. Then we introduce the fluxes: 


^^max 

no(A:) = / To{k)dk, 

Jk 

(4.9) 

r ^max 

Il\Y{k) = / Tw{k)dk^ 

Jk 

(4.10) 


and similarly for the perpendicular and parallel fluxes ns(A;_L) and ns(fc||). 

The isotropic energy fluxes for the slow modes and for the waves are shown in Fig. 12 
for the latest time {t = 30). We see that for all N/ f ratios, the energy fluxes for the slow 
modes are negative, and more or less independent of k, at scales larger than the forcing. 
Furthermore, the magnitude of this negative flux seems to be larger for relatively small 
values of N/f (the runs with N/f = 2,3,4 have about the same magnitude), while it 
is reduced when stratification prevails. On the contrary, the fluxes are positive at scales 
smaller than the forcing. In this range, the positive flux is an increasing function of N/f 
On the other hand, the energy fluxes for the wave modes do not change sign at the 
forcing scale: they are always positive. The absence of a range with constant small-scale 
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Figure 12. Isotropic fluxes of energy for the slow modes (top) and for the wave modes (bottom) 
for the series of runs at constant Froude number (left) and at constant Rossby number (right) 
with various N/f. Note the different scale (they differ by a factor 10) for Do and IIiv. 


fluxes can be attributed to the small size of the k > kp wave number range, a problem 


remedied when using higher resolutions and/or lower forcing wavenumber (Pouquet & 
Marino 2013 Marino et aZ.|2015), or a parametrization scheme (Pouquet et aZ.|20T3|. Note 


that unlike the traditional (total) energy fluxes, we need not have ns(0) = 0, because 
Eq and E\y are not independently conserved (only their sum is). This means that care 
should be taken when interpreting fluxes, since the value of the flux for one given mode 
(slow or wave) at k is not only the amount of energy transferred from wave numbers 
smaller than k to wave numbers larger than k, but also includes the contributions from 
the other mode. The total flux of energy from the slow modes to the waves flo—vU; which 
coincides with nw(0), or equivalently —no(0) (up to numerical errors), is shown as a 
function of N/f in Fig. 13 It is small when N/f is small (up to N/f = 4) and then 
increases almost linearly with N/f. For all values of N/f larger than unity, the energy 
flux for the wave modes (Fig. 12) is never significantly larger than Ho—w, which means 
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- 10 -^ 



Figure 13. Total energy flux from the slow modes to the wave modes as a function of N/f. 


that the transfer of wave energy is dominated by the leakage from slow modes. However, 
note that for the two cases N/f = 1/4 and N/f = 1/2, the flux at scales smaller than 
the forcing scale is much larger than Ho-j-vu, meaning that in those cases, there is a 
significant transfer of energy towards the small-scales from the waves themselves. This 
interpretation is consistent with the fact that in the N/f = 1/4 case the energy in the 
wave modes is significantly larger than in any other case, and these two cases are those 
for which the waves dominate at small scales (see Fig. [^, and for which the buoyancy 
Reynolds number is significantly larger. Similarly, the negative flux at scales larger than 
the forcing scales for the slow modes is on the order of no_>.vu when N/f > 7. This 
means that the negative flux alone does not allow to claim with certainty that we have 
an inverse cascade of slow modes. For instance, slow energy in the stratified case has a 
negative flux while we have seen above that it lacks some characteristics of the inverse 
cascade. What is more meaningful is the fact that in the range of scales between the 
forcing scale and the plateau of the flux at the largest scales, the flux decreases with k 
and goes below the level of Ho-j-vu, all the more so that N/f is small, consistently with 
the above analysis. For the runs with N/f < 4, the interpretation of the negative flux of 
slow energy is much easier: the magnitude of the more or less constant flux over scales 
larger than the forcing scales is much larger than no_>.vu, which means that it can safely 
be attributed to an inverse cascade of slow mode energy. 

In Fig. |14[ we show the perpendicular fluxes of energy in the slow modes and in the 
waves, no(fc_L) and 71w{ki_). For the slow modes, the fluxes are similar to the isotropic 
fluxes, with a notable difference: the negative fluxes at scales larger than the forcing 
scale are enhanced. In particular, now all the runs have a wide range of scales (larger 
than the forcing scale) for which the flux has larger magnitude than Ho-j-vu. That means 
that regardless of N/f, slow mode energy tends to be transferred towards large horizontal 
scales. This is known to be true even for purely stratified flows or flows strongly dominated 
by rotation, which tend to form uniform layers (Riley & belong 2000). As a matter of 
fact, runs dominated by stratification have a (negative) peak in the slow mode energy flux 
whose magnitude increases with N/ f. The fluxes indicate, however, that for stratification 
dominated runs (e.g. N/f = 20 or oo), this transfer of slow energy is arrested at some 
horizontal scales (like kinetic energy; see Marino et al. (2014)), while it continues all the 
way to the largest scales when stratification is weaker (e.g. N/f = 2,3,4 or even 7 and 
10). On the contrary, the energy flux in the wave modes are always smaller than Ho-j-vv, 
and contrary to the isotropic fluxes, they do not have any range of constant value. This 
indicates that the transfer of energy from the slow modes to the waves probably occurs 
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Figure 14. Perpendicular fluxes of energy for the slow modes (top) and for the wave modes 
(bottom) for the series of runs at constant Fronde number (left) and at constant Rossby number 
(right) with various N/f. 


at large horizontal scales (we cannot in principle draw definitive conclusions since there 
might be cancellations between energy transfer from slow to wave modes and transfer 
across scales of wave modes, but it does not seem very plausible that they are of different 
signs). Note that the forcing scale is difficult to identify from this figure alone. It should 
be emphasized that since we are forcing in a spherical shell in Fourier space, energy is 
actually injected in the full range G [0, kp] (the same goes for fey of course). Hence a 
caveat in the interpretation of anisotropic fluxes is that the forcing scale is not so well 
defined. However, in almost all cases except the above mentioned, the forcing scale can 
still be identified unambiguously in n 5 (A:_L) and ns(fc||), which gives us confidence that 
the forcing is not contaminating the results too much. Finally, the perpendicular fluxes 
of wave energy for the two cases N/f = 1/4 and N/f = 1/2 are very similar to the 
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Figure 15. Parallel fluxes of energy for the slow modes (top) and for the wave modes (bottom) 
for the series of runs at constant Fronde number (left) and at constant Rossby number (right) 
with various N/f. 


isotropic fluxes, which means that in those cases, the near-inertial waves tend to transfer 
their energy towards small horizontal scales. 

Now, we turn to the examination of the parallel flux of energy in the slow modes and 
in the wave modes (Fig. [T^. The most notable results are that runs with N/f >7 have 
positive fluxes of slow energy across the whole range of vertical scales, corresponding to 
the maintenance of horizontal layers, whose thickness decreases with N/f. This process 
is less clear when stratification is weaker, since the parallel flux then remains negative 
at the very large scales. For all the runs with N/f > 2, the parallel flux of wave energy 
is identical to the isotropic flux, indicating that vertical transfers dominate, although 
interpretation is still polluted by the exchange with slow modes. When N/f = 1/4 or 
N/f = 1/2, it is interesting to note that there is a peak at vertical scales larger than 
the forcing scale (this is another case where the isotropic forcing scale does not appear 
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clearly), although its magnitude is relatively small. This peak does not correspond to 
any identifiable feature in the isotropic fluxes. 


5. Discussion 

5.1. Summary of the results 

To sum up, § |4.1[ [4^ |4.3| and |4.4| have established that in all the runs shown here: 

• The total energy is dominated by the slow mode component, even when stratification 
is strong, although the fraction of energy in the waves increases when N/ f > 2 increases. 
This happens in spite of the fact that the forcing has roughly equal projections on the slow 
and wave modes. The increase in wave energy with N/ f is consistent with the observed 
total flux of energy from the slow modes to the wave modes, which also increases with 

N/f. 

• The slow mode component dominates at all scales, except in the two runs with higher 
buoyancy Reynolds number (those with N/f = 1/4 and 1/2), where they only dominate 
at large scale. This also holds when distinguishing horizontal and vertical scales. 

• When N/f is small enough (roughly N/f < 7, although the specific value may 
depend on the details of the setup, and, in particular, the buoyancy Reynolds number 
and the strength of the waves due to rotation and/or stratification), the slow mode 
component has all the characteristics of an inverse cascade of energy: linear growth of 
their global energy, self-similar growth of their energy at scales larger than the forcing 
scales, and negative and constant fluxes at scales larger than the forcing scale. 

• When TV// is higher, this inverse cascade of energy is arrested. When stratification 
prevails, slow modes dominate the total energy, and there is some transfer of these modes 
towards the large scales, but it does not take the form of an inverse cascade: it saturates 
with a flat spectrum at large isotropic and vertical scales, and there is no significant 
negative flux of energy for those scales. The spectrum is not flat at large horizontal 
scales, consistent with the presence of significant negative flux of slow mode energy, but 
it is still much shallower than in the presence of rotation. 

• As long as rotation is present, the fraction of energy in horizontally homogeneous 
modes (VSHF modes, satisfying k± = 0) does not vary much with N/f, and this fraction 
is always much lower than the energy in the waves, which is itself much lower than the 
energy in the slow modes. On the contrary, in the purely stratified case, the VSHF modes 
reach a level of energy comparable with the waves. Interestingly, this transition appears 
to be pretty sharp, contrary to all the other properties of the flows we have studied here 
(e.g. global partition of energy between waves and slow modes, growth rate of the slow 
mode energy and spectral indices), which seem to have a smoother transition as N/f 
increases. 


5.2. The role of slow modes, waves and VSHF modes for the direction of the energy 

cascade 


One aspect of the results outlined above (§ 5.11 which is of particular interest is that the 


discrepancy between the phenomenologies of rotating-stratified flows (inverse cascade of 
energy, when there is sufficient rotation), and purely stratified flows (direct cascade of 
energy) does not seem to be explained only by the different role of waves but also by the 
intrinsic difference in the behaviour of the slow modes in these two cases. 

A framework which is quite robust for understanding direct and inverse cascades is 
that of the inviscid invariants of the system. More precisely, the existence of a second 
quadratic, definite positive invariant, in addition to energy, may prevent the downscale 
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cascade of energy and lead to an inverse cascade. This is what happens in 2D turbu¬ 


lence (Kraichnan & Montgomery 19801. In stratified turbulence, with or without rota¬ 


tion, there is such a second inviscid invariant, though not quadratic in principle: potential 


enstrophy (see § 2.3). Although this invariant exists regardless of whether there is rota¬ 


tion or not, the phenomenology is not the same. A subtlety here is that the invariant 
has constant sign (it is always positive), but is not definite: some combinations of the 
fields make potential vorticity vanish. Such combinations correspond to the wave modes. 
Hence, a natural way of thinking is to say that if waves and slow modes did not exchange 
energy, the wave and slow mode energy would be separately conserved, and slow modes 
would have invariants similar to 2D turbulence, while waves would have invariants similar 
to 3D Homogenous Isotropic Turbulence. Therefore, in that case one should expect an 
inverse cascade of slow mode energy and a direct cascade of wave energy. 

However, in the purely stratified case, this reasoning does not hold since then, potential 
enstrophy is still degenerate even when we restrict it to the slow modes. Indeed, in the 
non-rotating case, only slow modes with fc_L 7 ^ 0 contribute to potential enstrophy (see 


§ 3.3). This peculiarity was used in a statistical mechanics argument (Herbert et al. 


2014) in order to explain why the two cases have opposite cascade directions in spite of 


having the same inviscid invariants. From this point of view, it could also be because 
of the presence of the vertically sheared horizontal flow modes that stratified turbulence 
does not have an inverse cascade. A similar situation occurs in homogeneous isotropic 
turbulence, where the second invariant, helicity, is not sign-definite in general (in our 
case the sign of potential enstrophy is constant but it can still vanish for some modes). 
When one decimates the system by keeping only helical modes of a given sign, helicity 
becomes sign-definite and an inverse cascade results, as observed numerically ([Biferale 


et al. 2012) and in agreement with statistical mechanics arguments (Zhu et al. 2014 
[Herbert 2014). But when one starts reintroducing in the dynamics a fraction of the 


modes which destroy sign-definiteness of the invariant, the inverse cascade collapses very 
soon ( Sahoo et aZ.||2015 ), which again seems in agreement with the statistical mechanics 
prediction (Herbert 2014) and the behavior of purely stratified flows described above. 


An alternative constraint which leads to the reversal of the energy flux is confinement: 
when the aspect ratio of the domain is small, enstrophy production is inhibited at scales 
larger than the layer thickness, even in the presence of stratification (Smith et aZ.||1996 


Celani et al. 2010 Sozza et al. 2015). To further characterize a potential mechanism 


involving the VSHF modes which prevents the inverse cascade in the purely stratified 
case, it would be interesting to look at triadic exchanges in more details. This is beyond 
the scope of this paper, and it is left for future studies. 

5.3. Variation with buoyancy Reynolds number 
As pointed out by several authors (Brethouwer et al. 2007 Ivey et ar||2008 [Bartello "fc 


Tobias 2013 Waite 2013), the buoyancy Reynolds number TZb plays an important role 
in stratified turbulence, and it is likely doing so as well when rotation is included. In 
particular, the regime of low buoyancy Reynolds number corresponds to a case where 
horizontal layers are viscously coupled. It is characterized by a vertical scale set by 
the viscous scale yJvLh/U rather than the buoyancy scale U/N (Riley & de Bruyn Kops 


2003 Godoy-Diana et a/.|2004 Brethouwer et al.\2007 Rorai et a7|2014 ), a small vertical 

Froude number (on the contrary, at high buoyancy Reynolds number, the vertical Froude 
number remains of order unity ( [Billant fc Chomaz|2001 )) and a potential enstrophy which 
is dominated by its quadratic component ( Waite|2013 ). Here we have explicitly verified 


the latter property (see § 2.3). Hence, the results presented here tend to show that in 
this regime, the system relaxes to a balanced state where the velocity field is almost 
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horizontal (the slow modes have vanishing vertical velocity), and slow mode energy is 
cascaded upscale provided rotation is strong enough. 

For higher TZb, data from the two weakly stratified runs point to the fact that waves 
now dominate the small-scales. This is bound to happen when the Ozmidov scale is 
resolved, since then isotropy is recovered and there is no longer any reason for the velocity 
field to be mostly horizontal: since the vertical velocity in the normal mode decomposition 
framework is entirely accounted for by the wave modes, their contribution has to be 
larger at these scales. Note that we have seen in § |2.3| that in this regime, the nonlinear 
components of potential vorticity are also stronger (see Fig. [^. In this high TZb regime, 
numerical simulations have shown that both an inverse and a direct cascade of energy 
were possible (Pouquet & Marino 2013 Marino et al. 2015[ ). From the point of view 
of the inviscid invariants of the system, this is a puzzling result since inverse cascades 


usually arise when a second invariant prevents the direct cascade (see § 5.2). The key 


to this phenomenon presumably lies in the fact that the energy breaks up into a wave 
component (not constrained by potential enstrophy conservation, as explained above) 
and a vortical component (constrained by potential enstrophy conservation), although it 
could be expected that in this regime, more energy is transferred from the slow modes 
to the waves than in the case investigated in the present paper. However, as soon as 
one resolves the inverse cascade, the Reynolds number, and therefore TZb are limited to 
moderate values since the forcing wavenumber is necessarily relatively large. Hence, such 
a study can only be done at high resolution. We leave this analysis for future work. 


6. Conclusion 

We have studied the partition of energy in terms of waves and vortices, defined through 
the normal modes of the linearized equations, in stratified flows with or without rotation. 
We have focused on the inverse cascade regime, which appears when rotation is sufficiently 
strong but not when stratification dominates. Using numerical simulations covering a 
wide range of N/ f ratios, we have characterized this transition in terms of the waves and 
vortices. 

Our main finding is that in all cases, the total energy is dominated by the vortices, 
and in most cases (i.e. when the buoyancy Reynolds number is of order one) this is true 
at all scales, vertical or horizontal alike. We have shown that in the range of parameters 
where the inverse cascade is present (i.e. when the ratio N/f is small enough), the 
vortices themselves exhibit all the characteristics of the inverse cascade: linear increase 
of the global energy in time, self-similar growth of the energy spectrum at scales larger 
than the forcing scale, and negative and (approximately) constant energy fluxes in the 
inertial range. When TV// becomes larger (here because rotation is decreased), these 
features disappear, in agreement with the weakening of the inverse cascade observed 
previously ( Marino et a/.||20l'3 2014). 

This is a puzzling result in that the inviscid invariants, which are traditionally reliable 
indicators of the direction of the energy cascade, are the same, regardless of the N/ f ratio, 
including in the purely stratified case. There are several elements which may contribute 
to the weakening, and ultimately, disappearance, of the inverse cascade of slow modes 
as N/f increases. On the one hand, although vortices dominate in all cases, the total 


flux of energy from slow modes to waves increases as N/f increases (see Fig. 13), and 
consequently the fraction of energy in the waves also increases (see Fig.[^. Besides, in the 
purely stratified case, there is an accumulation of energy in the VSHF modes, which do 


not contribute to quadratic potential enstrophy (see Fig. 11). Finally, as N/f increases. 


the approximation of potential enstrophy by its quadratic part, on which the inviscid 
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invariants argument is based, becomes less and less accurate (see Fig.[^. Future research 
would be needed to fully disentangle the relative role of these three effects. 
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